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Making sense of palaeoclimate sensitivity 
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Many palaeoclimate studies have quantified pre- anthropogenic climate change to calculate climate sensitivity (equi- 
librium temperature change in response to radiative forcing change) , but a lack of consistent methodologies produces a 
wide range of estimates and hinders comparability of results. Here we present a stricter approach, to improve 
intercomparison of palaeoclimate sensitivity estimates in a manner compatible with equilibrium projections for future 
climate change. Over the past 65 million years, this reveals a climate sensitivity (in K W -1 m 2 ) of 0.3-1. 9 or 0.6-1.3 at 95% or 
68% probability, respectively. The latter implies a warming of 2.2-4. 8 K per doubling of atmospheric C0 2 , which agrees 
with IPCC estimates. 


C haracterizing the complex responses of climate to changes in the 
radiation budget requires the definition of climate sensitivity: 
this is the global equilibrium surface temperature response to 
changes in radiative forcing (an alteration to the balance of incoming and 
outgoing energy in the Earth-atmosphere system) caused by a doubling 
of atmospheric C0 2 concentrations. Despite progress in modelling and 
data acquisition, uncertainties remain regarding the exact value of cli- 
mate sensitivity and its potential variability through time. The range of 
climate sensitivities in climate models used for Intergovernmental Panel 
for Climate Change Assessment Report 4 (IPCC-AR4) is 2.1-4.4K per 
C0 2 doubling 1 , or a warming of 0.6- 1. 2 K per W m -2 of forcing. 
Observational studies have not narrowed this range, and the upper limit 
is particularly difficult to estimate 2 . 

Large palaeoclimate changes can be used to estimate climate sensi- 
tivity on centennial to multi-millennial timescales, when estimates of 
both global mean temperature and radiative perturbations linked with 
slow components of the climate system (for example, carbon cycle, land 
ice) are available (Fig. 1). Here we evaluate published estimates of climate 
sensitivity from a variety of geological episodes, but find that intercom- 
parison is hindered by differences in the definition of climate sensitivity 
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Figure 1 | Typical timescales of different feedbacks relevant to equilibrium 
climate sensitivity, as discussed in this work. Modified and extended from 
previous work 98 . Ocean timescales were extended to multi-millennial timescales 99 . 


between studies (Table 1). There is a clear need for consistent definition 
of which processes are included and excluded in the estimated sensitivity, 
like the need for strict taxonomy in biology. The definition must agree 
as closely as possible with that used in modelling studies of past and 
future climate, while remaining sufficiently pragmatic (operational) to 
be applicable to studies of different climate states in the geological past. 

Here we propose a consistent operational definition for palaeoclimate 
sensitivity and illustrate how a tighter definition narrows the range of 
reported estimates. Consistent intercomparison is crucial to detect sys- 
tematic differences in sensitivity values — for example, due to changing 
continental configurations, different climate background states, and the 
types of radiative perturbations considered. These differences may then 
be evaluated in terms of additional controls on climate sensitivity, such 
as those arising from plate tectonics, weathering cycles, changes in 
ocean circulation, non-C0 2 greenhouse gases (GHGs), enhanced water- 
vapour and cloud feedbacks under warm climate states. Palaeoclimate 
data allow such investigations across geological episodes with very dif- 
ferent climates, both warmer and colder than today. Clarifying the 
dependence of feedbacks, and therefore climate sensitivity, on the back- 
ground climate state is a top priority, because it is central to the utility of 
past climate sensitivity estimates in assessing the credibility of future 
climate projections 1,3 . 

Quantifying climate sensitivity 

‘Equilibrium climate sensitivity’ is classically defined as the simulated 
global mean surface air temperature increase (AT, in K) in response to a 
doubling of atmospheric C0 2 , starting from pre- industrial conditions 
(which corresponds to a radiative perturbation, A R, of 3.7 W m -2 ; refs 1, 
3). We introduce the more general definition of the ‘climate sensitivity 
parameter’ as the mean surface temperature response to any radiative 
perturbation ( S = AT/AR ; where AT and A R are centennial to multi- 
millennial averages), which facilitates comparisons between studies 
from different time- slices in Earth history. For brevity, we refer to S as 
‘climate sensitivity’. In the definition of S, an initial perturbation AR 0 
leads to a temperature response AT 0 following the Stefan-Boltzmann 
law, which is the temperature-dependent blackbody radiation response. 
This is often referred to as the Planck response 4 , with a value S 0 of about 
0.3KW” 1 m 2 for the present-day climate 5,6 . The radiative perturbation 
of the climate system is increased (weakened) by various positive (nega- 
tive) feedback processes, which operate at a range of different timescales 
(Fig. 1). Because the net effect of positive feedbacks is found to be greater 
than that of negative feedbacks, the end result is an increased climate 
sensitivity relative to the Planck response 4 . 


*Lists of participants and their affiliations appear at the end of the paper. 
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Table 1 | Summary of key studies. 


Label in Fig. 3 

Source 

Time window 

Explicitly considered 

Temperature data used 

S and la bounds 

Notes 




forcings 


(KW _1 m 2 ) 


1 

Ref. 2 

LGM 

Various 

Various 

0.81 ±0.27 

LGM compilation based 






(data); 

0.81 +§ 2 7 (models) 

on ref. 15 

2 

Ref. 6 

LGM 

GHG(C0 2 ,CH 4 , N 2 0), 

ATg| 0 bai = —5.8 ± 1.4 K; 

0 79 + °- 33 
u - ' -0.23 

Scaling factor (0.85) for 




LI, AE, VG 

GLAMAP extrapolated 

smaller S at LGM compared 





with model 82 


to 2 x C0 2 (refs 12, 16) 

3 

Ref. 86 

LGM 

GHG(C0 2 ,CH 4 ), 

CLIMAP and AT aa&g | d 

0.80 ±0.14 

Value after authors’ suggested 




LI, AE, VG 


correction of CLIMAP 
temperatures 


4 

Ref. 79 

LGM 

GHG (C0 2 ,CH 4 ), 

MARGO 81 SST based 

0.62+° ; ° 8 

Model-based global estimate 




LI, AE, VG 

A T= -3.0+J; 8 K 



5 

Ref. 76 

GC 

GHG (C0 2 ,CH 4 ) 

A7"tr 0 p 

1.1 ±0.05 

Author’s linear regression case. 
Value based on single-site tropical 
SST, and representation of global 
changes will be more uncertain 

6 

Ref. 74 

GC 

GHG (C0 2 ,CH 4 ), 

A T aa (with 1.5 x polar 

0.88 ±0.13 

Author used a single value for 




LI, AE 

amplification) 


polar amplification. If 2 x were 
used 52 , then the central estimate 
is closer to 0.7 

7 

Ref. 52 

GC 

GHG (C0 2 ,CH 4 , 

A T aa (with 2 x polar 

0.75 ±0.13 

Authors used a single value for 




N 2 0), LI 

amplification) 


polar amplification. If 1.5 x were 
used 74 , then the central estimate 
becomes 1.0 

8 

Ref. 52 

GC 

GHG 

A T aa (with 2 x polar 

1.5 ±0.25 

Authors used a single value for 




(C0 2 , CH 4 , N 2 0) 

amplification) 


polar amplification. If 1.5 x were 
used 74 , then the central estimate 
becomes 2.0 

23-32 

This work, 

GC (<800 kyr ago) 

GHG (C0 2 , CH 4 , N 2 0), 

AT nh = model-based 

0.66 ± 0.22 to 

This covers the range of 


based on ref. 6 


LI, AE, VG 

deconvolution of benthic 
8 1S 0 (ref. 51), scaled to 
global AT using a NH 
polar amplification on 
land of 2.75 ±0.25 

2.26 ±0.78 

S[GHG,X] S iven in Table 2 



9 

Ref. 85 

GC 

GHG (C0 2 ,CH 4 , 

A T aa (with 2 x polar 

0.75 ±0.13 

Authors used a single value for 




N 2 0), LI 

amplification) and 


polar amplification. If 1.5 x were 





1.5 x A 7 ds 


used 73 , then the central estimate 
becomes 1.0 

10 

Ref. 39 

GC 

GHG(C0 2 ,CH 4 , 

36-record global SST 

0.85+gJ 5 

Polar amplification diagnosed, 




N 2 0), LI, AE 

synthesis along with 

not imposed. Estimates made 





AT aa&g | d . 


both in a spatially explicit sense 
and as direct global means 

11 

Ref. 39 

GC 

GHG(C0 2 ,CH 4 , 

36-record global SST 

1.05 ±0.25 

As above 




N 2 0), LI 

synthesis along with 
AT aa&g | d . 



12 

Ref. 87 

Early to Middle 

C0 2 , ESS 

Using model-based 

1.92 ±0.14 to 

Forcing in ref. 44; temperature 



Pliocene (4.2-3.3 


AT for Middle and 

2.35 ±0.18 

in ref. 87. Both derived in global 



Myr ago) 


Early Pliocene of 

(3.3 Myr ago); 

sense from model experiments 





2.4-2.9 °C and 4 °C. 

2.60 ±0.19 






AC0 2 alkenone 

(4.2 Myr ago) 


13 

Ref. 61 

Miocene optimum 

Slow feedbacks 

Deconvolution of 

0.78 ± 10% 

f= 0.71, £ = 5.35, y= 1.3. Details 



to present day 


benthic S 18 0 (ref. 63) 


in Supplementary Information 

14 

This work 

Eocene-Oligocene 

C0 2 . ESS (in the 

Model-based AT, 

1 70 + 0.9 
1 - / ^-0.54 

Details in Supplementary 


(compilation) 

transition (—34 Myr 

sense of ref. 44) 

with range of C0 2 

Information 



ago) 


values 



15 

This work 

Late Eocene 

C0 2 . ESS (in the 

Model-based AT, 

1 R? + 0.26 

Details in Supplementary 


(compilation) 

versus present 

sense of ref. 44) 

with range of C0 2 values 

Information 

16 

Ref. 78 

Middle Eocene 

C0 2 . Ice-free world. 

A7 ds (2 records) and 

0.95 ± 0.3 

500 kyr timescale. A7 ds = A7 mg . 



Climatic Optimum 

Event study (not 

A 7 mg (7 records). 


Temperatures from subtropics to 



(-40 Myr ago) 

affected by plate 

AC0 2 from alkenones 


high latitudes; no tropical data. 




tectonics and 



Hence biased to high-latitude 




evolution effects) 



sensitivity 

17 

Ref. 78 

Mid to Late Eocene 

C0 2 . Largely ice-free 

A 7 ds (ref. 71) and A T mg . 

0.95 ± 0.3 

Multi-million-year timescale. 



transition (41-35 

world. Event study 

AC0 2 = difference mid 


Adding uncertainty of ± 1 °C 



Myr ago) 

(not affected by 

Eocene alkenone and 


to AT would enhance 1 0 limits 




plate tectonics and 
evolution effects) 

late Eocene S n B 


to ±0.45 KW -1 m 2 

18 

Ref. 88 

Early Eocene 

C0 2 . Ice-free world. 

AT mg (refs 89-91). AC0 2 

0.65 ± 0.25 

Central value recalculated in 



(-55-50 Myr ago) 

(potential influences 

based on modelling 91 


ref. 94. 




of plate tectonics and 

marine organic carbon 


Note ref. 89 underestimated 




biological evolution 

isotope fractionation 92 


tropical SST 




not considered) 

and soil nodules 93 



19 

This work 

PETM 

C0 2 . Ice-free world. 

A 7 ds (>6 records) and 

1.0-1. 8 

Details in Supplementary 


(compilation) 

(-56 Myr ago) 

Event study (not 

A T mg (>1 1 records; 


Information. Assumes all warming 




affected by plate 

equatorial to polar). 


due to C input, and range of 




tectonics and 

AC0 2 based on deep ocean 

background C0 2 and C-injection 




evolution effects) 

carbonate chemistry 72 ’ 95 


scenarios. A7 ds = A 7 mg . Total 
range of S is 0.7-2.2 K W^ 1 m 2 . 
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Table 1 | Continued 


Label in Fig. 3 Source 

Time window 

Explicitly considered 
forcings 

Temperature data used 

S and ler bounds 
(KW -1 m 2 ) 

Notes 

20 

Ref. 96 

Cretaceous and 
early Palaeogene 

C0 2 . Largely ice-free 
world, (potential 
influences of plate 
tectonics and biological 
evolution not 
considered) 


1 

Recalculated in ref. 94. No 
uncertainty range was reported, 
nor salient details for assessment. 
Figure 3b, c assumes ±25% 

21 

Ref. 94 

Cretaceous and 
early Palaeogene 

C0 2 . Largely ice-free 
world. ESS in the 
sense of ref. 44 

AT after refs 52, 71. 
AC0 2 based on ref. 60. 

>0.8 

No uncertainty range reported. 
This is a lower bound estimate 
only 

22 

Ref. 97 

Phanerozoic 

C0 2 . Ice-free situation. 
(Potential influences 
of plate tectonics and 
biological evolution 
not considered). 

A7" mg AC0 2 based 
on GEOCARBSULF 

0.8-1.08 

Model-based with extensive 
uncertainty analysis 


These studies have empirically determined S for the Pleistocene and some deep-time periods from comparison between data-derived time series for temperature and for radiative change. Comparison of results 
between studies is greatly hindered by the different ‘versions’ of S used, as related to different notions of which processes should be explicitly accounted for, and by the different approaches taken to approximate 
global mean surface temperature. All uncertainties areas originally reported, but shown here at the level equivalent to la, estimated where necessary by dividing total range values by a factor of 2. All values for Sare 
reported in KW -1 m 2 , where necessary after transformation using 3.7 W rrf 2 per doubling of CO 2 , bearing in mind the caveats for this at high CO 2 concentrations as elaborated in the main text. GC, glacial cycles; 
LGM, Last Glacial Maximum; PETM, Palaeocene/Eocene thermal maximum; SST, sea surface temperature. See main text for details of forcings. Subscripts: aa, Antarctica; gld, Greenland; trop, tropical; ds, deep 
sea; global, global mean; mg, Mg/Ca; NH, Northern Hemisphere. 


We emphasize that all feedbacks, and thus the calculated climate 
sensitivity, may depend in a — largely unknown — nonlinear manner 
on the state of the system before perturbation (the ‘background climate 
state’) and on the type of forcing 7-15 . The relationship of S with back- 
ground climate state differs among climate models 12,16-18 . A suggestion 
of state dependence is also found in a data comparison (Table 2) 6 , where 
climate sensitivity for the past 800,000 years (800 kyr) shows substantial 
fluctuations through time (Fig. 2). In contrast, its values for the Last 
Glacial Maximum (LGM) alone occupy only the lower half of that range 
(Fig. 2). That evaluation also suggests that the relationship of S with the 
general climate state may not be simple. 

‘Fast' versus ‘slow' processes 

Climate sensitivity depends on processes that operate on many different 
timescales, from seconds to millions of years, due to both direct response 
to external radiative forcing, and internal feedback processes (Fig. 1). 
Hence, the timescale over which climate sensitivity is considered is 
critical. An operationally pragmatic decision is needed to categorize a 
process as ‘slow’ or ‘fast’, depending on the timescale of interest, the 
resolution of the (palaeo-) records considered and the character of 
changes therein 19 . If a process results in temperature changes that reach 
steady state slower than the timescale of the underlying radiative per- 
turbation, then it is considered ‘slow’; if it is faster or coincident, then it is 
‘fast’. Present-day atmospheric GHG concentrations and the radiative 
perturbation due to anthropogenic emissions increase much faster than 
observed for any natural process within the Cenozoic era 20-22 . 


For the present, the relevant timescale for distinguishing between fast 
and slow processes can be taken as 100 yr (ref. 23). Ocean heat uptake 
plays out over multiple centuries. Combined with further ‘slow’ pro- 
cesses, it causes climate change over the next few decades to centuries to 
be dominated by the so-called ‘transient climate response’ to radiative 
changes that result from changing GHG concentrations and aerosols 5,19,24 . 
After about 100 yr, this transient climate response is thought to amount to 
roughly two-thirds of the equilibrium (see below) climate sensitivity 5,25 . 
Climate models account for the fast feedbacks from changes in water- 
vapour content, lapse rate, cloud cover, snow and sea-ice albedo 26 , and 
the resulting response is often referred to as the ‘fast- feedback’ or 
‘Charney’ sensitivity 23 . To approximate the ‘equilibrium’ value of that 
climate sensitivity, accounting for ocean heat uptake and further slow 
processes, models might be run over centuries with all the associated 
computational difficulties 27-30 , or alternative approaches may be used 
that exploit the energy balance of the system for known forcing or 
extrapolation to equilibrium 31 . 

In palaeoclimate studies, an operational distinction has emerged to 
distinguish ‘fast’ and ‘slow’ processes relative to the timescales of tempera- 
ture responses measured in palaeodata, where ‘fast’ is taken to apply to 
processes up to centennial scales, and ‘slow’ to processes with timescales 
close to millennial or longer. Thus, changes in natural GHG concentrations 
are dominated by ‘slow’ feedbacks related to global biogeochemical cycles 
(Fig. 1). Similarly slow are the radiative influences of albedo feedbacks that 
are dominated by centennial-scale or longer changes in global vegetation 
cover and global ice area/ volume (continental ice sheets) (Fig. 1). 


Table 2 | Common permutations of S that may be encountered in palaeostudies 


Label in 

S definition 

Explicitly considered 

Period in which it is practical 

S± la for 800 kyr 

S± la 

S 

Fig. 3 


radiative perturbation 

to use the definition 

(KW -1 m 2 ) 

for LGM (KW -1 m 2 ) 

for Pliocene (KW -1 m 2 ) 

23 

CO 

o 

o 

ro 

Aft[C02] 

All (especially pre-35 Myr ago 
when LI ~ 0) 

3.08 ± 0.96 

2.63 ±0.57 

1.2 

24 

S[C02, LI] 

A^[C02, LI] 

<35 Myr ago 

1.07 ±0.40 

0.95 ±0.22 

0.97 

25 

S[C02, LI, VG] 

Aft[C02, LI, VG] 

<35 Myr ago 

0.86 ±0.27 

0.80 ±0.19 

0.82 

26 

S[C02, LI, AE] 

Aft[C02, LI, AE] 

<35 Myr ago, but mainly <800 kyr ago 

0.90 ± 0.42 

0.72 ±0.18 


27 

S[C02, LI, AE, VG] 

Aft[C02, LI, AE, VG] 

<35 Myr ago, but mainly <800 kyr ago 

0.75 ± 0.29 

0.63 ±0.15 


28 

S[GHG] 

A/?[GHG] 

<800 kyr ago 

2.32 ± 0.76 

1.97 ±0.41 


29 

S[GHG, LI] 

A%3HG, LI] 

<800 kyr ago 

0.96 ± 0.36 

0.85 ±0.19 


30 

S[GHG, LI, VG] 

AR*[GHG, LI, VG] 

<800 kyr ago 

0.78 ± 0.23 

0.73 ±0.16 


31 

S[GHG, LI, AE] 

A%3HG, LI, AE] 

<800 kyr ago 

0.82 ± 0.36 

0.66 ±0.16 


32 

S[GHG, LI, AE, VG] 

A%3HG, LI, AE, VG] 

<800 kyr ago 

0.68 ± 0.24 

0.58 ±0.14 



S (second column) is presented with a subscript that identifies the explicitly considered radiative perturbations A R (third column, same subscripts as for S); all other processes are implicitly resolved as feedbacks 
within S. The period in which the various definitions of S are practical is determined by the availability of data for the explicitly considered processes. Subscript C02 indicates the radiative impact of atmospheric 
C0 2 concentration changes; LI represents the radiative impact of global land ice-volume changes; VG stands for the radiative impact of global vegetation cover changes; AE indicates the radiative impact of aerosol 
changes; GHG stands for the impact of changes in all non-water natural greenhouse gases (notably C0 2 , CH 4 and N 2 0). Columns 5 and 6 give calculated values for all suggested permutations of S for the past 
800 kyr or the LGM, respectively, based on a previous data compilation 6 . Mean values of all S [X ]forthe LGM are about 13% smaller than for the whole 800 kyr, but Newell within the given uncertainties. This offset 
illustrates the state-dependence of S (see Supplementary Information). Column 7 gives examples for the Pliocene 13,44 ; Fig. 3b, c assumes ±25% uncertainty in these. In these values the effects of orographic 
changes have been taken into account (see Supplementary Information section B2). 
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Figure 2 | Illustration of variability of climate sensitivity using a calculation 
of S[co2,li]> as defined in this work, for the past 800 kyr. a, Changes in global 
temperature, b, Changes in radiative forcing due to changes in C0 2 and surface 
albedo due to land ice. c, Calculated S[co2,li]> which is only considered robust 
and calculated when AT < — 1.5 K and AT[ C o 2 ,li] < — 0.5 W m -2 , as indicated 
by the dotted red lines in a and b. In c, mean of ± a 0 (dashed black lines 
indicate a 0 , the uncertainty of averaging) and 100-kyr running mean (blue line) 
are shown. Magenta marker in c denotes S { ± o x for the LGM only (23-19 kyr 
ago) (or is the square root of the sum of squares of individual uncertainties 
connected with different processes contributing to S z ). The grey areas in 
a-c denote o x (standard deviation) uncertainties of S 2 for single points in time 
(points themselves are omitted for clarity). Details of data and the definition of 
the calculated uncertainties presented in this figure are available in 
Supplementary Information. In a and b, the dashed black lines indicate the 
preindustrial reference case (AT = 0 K, AT[ C o2,li] = 0 W m -2 ). 

Other processes clearly have both fast and slow components. For 
example, palaeorecords of atmospheric dust deposition imply important 
aerosol variations on decadal to astronomical (orbital) timescales 32-36 , 
reflecting both slow controlling processes related to ice-volume and 
land- surface changes, and fast processes related to changes in atmo- 
spheric circulation. A further complication arises from the lack of ade- 
quate global atmospheric dust data for any geological episode except the 
LGM 37,38 , even though that is essential because the spatial distribu- 
tion of dust in the atmosphere tends to be inhomogeneous and because 
temporal variations in some locations tend to take place over several orders 
of magnitude 32-36 . Moreover, palaeoclimate models generally struggle to 
account for aerosols, with experiments neither prescribing nor implicitly 
resolving aerosol influences. So far, understanding of aerosol/dust feed- 
backs remains weak and in need of improvements in both data coverage 


and process modelling, especially because dust forcing may account for 
some 20% of the glacial-interglacial change in the radiative budget 6,39 . 

So for comparison of results between studies, it is most effective to 
consider only the classical ‘Charney’ water- vapour, cloud, lapse rate, and 
snow and sea-ice feedbacks 23 as ‘fast’, and all other feedbacks as ‘slow’. In 
addition, results from palaeoclimate sensitivity studies generally do not 
address the transient climate response that dominates present-day 
changes, but capture a more complete longer-term system response 
comparable with equilibrium climate sensitivity in climate models. 

Forcing and slow feedbacks 

The external drivers of past natural climate changes mainly resulted 
from changes in solar luminosity over time 40 , from temporal and spatial 
variations in insolation due to changes in astronomical parameters 41-43 , 
from changes in continental configurations 14,44 , and from geological 
processes that directly affect the carbon cycle (for example, volcanic 
outgassing). However, the complete Earth system response to such forcings 
as recorded by palaeodata cannot be immediately deduced from the 
(equilibrium) ‘fast feedback’ sensitivity of climate models, because of 
the inclusion of slow feedback contributions. When estimating climate 
sensitivity from palaeodata, agreement is therefore needed about which 
of the slower feedback processes are viewed as feedbacks (implicitly 
accounted for in S), and which are best considered as radiative forcings 
(explicitly accounted for in A R). 

We employ an operational distinction 31,45 in which a process is con- 
sidered as a radiative forcing if its radiative influence is not changing 
with temperature on the timescale considered, and as a feedback if its 
impact on the radiation balance is affected by temperature changes on 
that timescale. For example, the radiative impacts of GHG changes over 
the past 800 kyr may be derived from concentration measurements of 
C0 2 , CH 4 and N 2 0 in ice cores 46-48 , and the radiative impacts of land-ice 
albedo changes may be calculated from continental ice-sheet estimates, 
mainly based on sea-level records 49-51 . Thus, the impacts of these slow 
feedbacks can be explicitly accounted for before climate sensitivity is 
calculated. This leaves only fast feedbacks to be considered implicitly in 
the calculated climate sensitivity, which so approximates the (equilib- 
rium) ‘Charney’ sensitivity from modelling studies 6,39,52 . 

Operational challenges 

All palaeoclimate sensitivity studies are affected by limitations of data 
availability. Below we discuss such limitations to reconstructions of 
forcings and feedbacks, and of global surface temperature responses. 
First, however, we re-iterate a critical caveat, namely that the climate 
response depends to some degree on the type of forcing (for example, 
shortwave versus longwave, surface versus top -of- atmosphere, and local 
versus global). The various radiative forcings with similar absolute mag- 
nitudes have different spatial distributions and physics, so that the con- 
cept of global mean radiative forcing is a simplification that introduces 
some (difficult to quantify) uncertainty. 

Astronomical (orbital) forcing is a key driver of climate change. In 
global annual mean calculations of radiative change, astronomical forcing 
is very small and often ignored 39,52 . Although this obscures its importance, 
mainly concerning seasonal changes in the spatial distribution of inso- 
lation over the planet 41,42,53-55 , we propose that the contribution of the 
astronomical forcing to A R may be neglected initially. When other com- 
ponents of the system respond to the seasonal aspects of forcing, such as 
Quaternary ice-sheet variations, these may be accounted for as forcings 
themselves. 

GHG concentrations from ice cores are not available for times before 
800 kyr ago, when C0 2 levels instead have to be estimated from indirect 
methods. These employ physico-chemical or biological processes that 
depend on C0 2 concentrations, such as the abundance of stomata on 
fossil leaves 56 , fractionation of stable carbon isotopes by marine phyto- 
plankton 57 , boron speciation and isotopic fractionation in sea water as a 
function of pH and preserved in biogenic calcite 58 , and the stability fields 
of minerals precipitated from waters in contact with the atmosphere 59 . 
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Considerable uncertainties remain in such reconstructions, but improve- 
ments are continually made to the methods, their temporal coverage and 
their mutual consistency 60 . Recent work has synthesized a high-resolu- 
tion C0 2 record for the past 20 million years (Myr; ref. 61), but new data 
and updated syntheses remain essential, particularly for warmer climate 
states. Also, proxies are needed for reconstruction of CH 4 and N 2 0 con- 
centrations in periods pre-dating the ice-core records 62 . 

Regarding the assessment of land- ice albedo changes, good methods 
exist for the generation of continuous centennial- to millennial- scale 
sea-level (ice-volume) records over the past 500 kyr (refs 49-51), but 
such detailed information remains scarce for older periods. A model- 
based deconvolution of deep-sea stable oxygen isotope records into 
their ice-volume and deep-sea temperature components 51 was recently 
extended to 35 Myr ago 63 , but urgently requires independent validation, 
especially to address uncertainties about the volume-to-area relation- 
ships that would be different for incipient ice sheets than for mature 
ice sheets 64,65 . Before 35 Myr ago, there is thought to have been (virtually) 
no significant land-ice volume 66 , but this does not exclude the potential 
existence of major semi-permanent snow/ice- fields 67,68 , and there remain 
questions whether these would constitute ‘fast’ (snow) or ‘slow’ (land-ice) 
feedbacks. The contribution of the sea-ice albedo feedback also remains 
uncertain, with little quantitative information beyond the LGM. 

Similar examples of uncertainties and limited data availability could 
be listed for all feedbacks. However, a ‘deep-time’ (before 1 Myr ago) 
geological perspective must be maintained because it offers access to the 
nearest natural approximations of the current rate and magnitude of 
GHG emissions 69,70 , and because only ancient records provide insight 
into climate states globally warmer than the present. Given that no past 
perturbation will ever present a perfect analogue for the continuing 
anthropogenic perturbation, it may be more useful to consider past 
warm climate states as test-beds for evaluating processes and responses, 
and for challenging/ validating model simulations of those past climate 
states. Such data-model comparisons will drive model skill and under- 
standing of processes, improving confidence in future multi-century 
projections. For such an approach, palaeostudies may minimize the 
impacts of very long-term influences on climate sensitivity (for example, 
due to changes in orography, or biological evolution of vegetation) 
through a focus on highly resolved documentation of specific perturba- 
tions that are superimposed upon different long-term background climate 
states. An example is the pronounced transient global warming and 
carbon-cycle perturbation during the Palaeocene/Eocene thermal maxi- 
mum (PETM) anomaly 71,72 , which punctuated an already warm climate 
state 73 . Note that deep-time case studies need to consider one further 
complication, namely that the radiative forcing per C0 2 doubling may 
be about 3.7Wm" 2 when starting from pre-industrial concentrations, 
but increases at higher C0 2 levels 11 . Data-led studies may help with a 
first-order documentation of this dependence. Calculation of S from C0 2 
and temperature measurements using a constant 3.7Wm" 2 per C0 2 
doubling would (knowingly) overestimate S for high-C0 2 episodes. 
The difference with other, identically defined, S values for different cli- 
mate background states may then be used to assess any deviation from 
3.7Wm' 2 per C0 2 doubling. 

Regarding the reconstruction of past global surface temperature res- 
ponses (that is, AT in equation (1) below), again much remains to be 
improved. Most work to date (see Table 1) relies on one or more of the 
following: polar temperature variations from Antarctic ice cores (since 
800 kyr ago) with a multiplicative correction for ‘polar amplification’ 
(usually estimated at 1. 5-2.0; refs 74, 75); deep-sea temperature varia- 
tions from marine sediment-core data with a correction for the ratio 
between global surface temperature and deep-sea temperature changes 
(often estimated at 1.5); single-site sea surface temperature (SST) records 
from marine sediment cores; or compilations of SST data of varying 
geographic coverage from marine sediment cores 6,39,52,76-78 . So far, few 
studies have included terrestrial temperature proxy records other than 
those from ice cores 79 , yet better control on land- surface data is crucial 
because of seasonal and land- sea contrasts. Continued development is 


needed of independently validated (multi-proxy) and spatially represent- 
ative (global) data sets of high temporal resolution relative to the climate 
perturbations studied. 

Uncertainties in individual reconstructions of temperature change 
may in exceptional cases be reported to ± 0.5 K, but more comprehensive 
uncertainty assessments normally find them to be larger 80,81 . Compilation 
of such records to determine changes in global mean surface tempera- 
ture involves the propagation of further assumptions/uncertainties, for 
example due to interpolation from limited spatial coverage, and the end 
result is unlikely to be constrained within narrower limits than ± 1 K even 
for well- studied intervals. Finally, comparisons between independent 
reconstructions for the same episode reveal ‘hidden’ uncertainties due 
to differences between each study’s methodological choices, uncertainty 
determination, and data- quality criteria, which are hard to quantify and 
often poorly elucidated. Take the LGM for example, which for tempera- 
ture is among the best-studied intervals. The MARGO compilation 81 
inferred a global SST reduction of -1.9 ± 1.8 K relative to the present. 
Another spatially explicit study 79 used that range to infer a global mean 
surface air temperature anomaly of — 3 1 J; 2 K. The latter contrasts with a 
previous estimate of — 5.8 ± 1.4 K (ref. 82), which is consistent with 
tropical (30° S to 30° N) SST anomalies of —2.7 ± 1.4 K (ref. 83). 
However, that tropical range itself is also contested; the MARGO 81 
study suggested such cooling in the Atlantic Ocean, but less in the 
tropics of the Indian and Pacific Oceans (giving a global tropical 
cooling of only — 1.7 ± 1.0 K). Clearly, even a well-studied interval 
gives rise to a range of estimates for temperature, and therefore for 
climate sensitivity. 

It is evident that progress in quantifying palaeoclimate sensitivity 
will not only rely on a common concept and terminology that allows 
like-for-like comparisons (see below); it will also rely on an objective, 
transparent and hence reproducible discussion in each study of the 
assumptions and uncertainties that affect the values determined for 
change in both temperature and radiative forcing. 

A way forward 

Here we propose a new terminology to help palaeoclimate sensitivity stud- 
ies adopt common concepts and approaches, and thus improve the poten- 
tial for like-for-like comparisons between studies. First we outline how our 
concept of ‘equilibrium’ S for palaeo- studies relates to ‘equilibrium’ S for 
modern studies. Then, we present a notation system that is primarily of 
value to data-based palaeo studies to clarify which slow feedbacks are 
explicitly accounted for. We finish with an application of the new frame- 
work, calculating climate sensitivity from a representative selection of 
palaeoclimate sensitivity estimates over the past 65 Myr, with a fair balance 
of climates warmer than the present to those colder than the present. 

When the AT response to an applied GHG radiative forcing A R is 
small relative to ‘pre-perturbation’ reference temperature, the ‘equilib- 
rium’ climate sensitivity S a (where a indicates actuo , for present-day) 
takes the form (see, for example, refs 4, 84): 

^4— (i) 


A R 


h + 


Here 2 P is the Planck feedback parameter ( — 3.2 W m -2 K -1 ) and 2 f z - 
(in W m -2 K -1 ) represents the feedback parameters of any number 
(N) of fast (f ) feedbacks. We define feedback parameters in the form 
2 f z = ARJ AT. S a is the ‘Charney’ sensitivity calculated by most cli- 
mate models in ‘2 X C0 2 ’ equilibrium simulations, with a range of 
0.6- 1.2 K W -1 m 2 in IPCC-AR4. However, the Earth system in reality 
responds to a perturbation according to an equilibrium climate sen- 
sitivity parameter Sf* (where p indicates palaeo ), but the timescales to 
reach this equilibrium are long, so that the forcing normally changes 
before equilibrium is reached. To obtain S a from palaeoclimate sens- 
itivity S p , a correction is therefore needed for the slow feedback influ- 
ences. Using X s j to represent any number (M) of slow (s) feedbacks, we 
derive the general expression (see Supplementary Information): 
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This approach is contingent on the above-mentioned caveats of state- 
dependence, linearization (small AT), changes in slow feedbacks, and 
transient effects, where the last is relevant only in records of excep- 
tionally high temporal resolution. Knowledge of slow (2 s ) and fast (2 f ) 
feedbacks can be combined into a factor F = 2 s /(2 f + 2 s ) that may 
then be used to back- calculate fast feedbacks out of palaeoclimate 
sensitivity S p . 

A recent study 44 defined the term ‘Earth system sensitivity’ (ESS) to 
represent the long-term climate response of Earth’s climate system to a 
given C0 2 forcing, including both fast and slow processes. In our nota- 
tion, ESS = AT 2 xco 2 where AT 2X co 2 is the forcing due to a C0 2 - 
doubling (3.7 W m -2 ). 

Here we introduce a more explicit notation regarding what was (not) 
included in the climate sensitivity diagnosis. It is the ‘specific climate 
sensitivity’ S[ A , B ...]> expressed in KW” 1 m 2 , where slow feedback pro- 
cesses A, B, and so on, are explicitly accounted for (that is, included in 
the forcing term, AT[ A ,b...]). We use ‘LI’ to denote albedo forcing due to 
land-ice volume/area changes, ‘VG’ for vegetation- albedo forcing, ‘AE’ 
for aerosol forcing and ‘C02’ for atmospheric C0 2 forcing (see also 
Table 1). This approach requires from the outset that a comprehensive 
view is taken of the various causes of change in the radiative balance. 


The most practical version of S to be estimated from palaeodata 
is S [C0 2 ,li]> because S [C0 2 ,li] = S[co 2 ] during times (pre-35Myr ago) 
without ice volume, and because global vegetation cover changes, atmo- 
spheric dust fluctuations, and both CH 4 and N 2 0 fluctuations (the two 
important non-C0 2 GHGs) generally remain poorly constrained by 
proxy data. Common reporting of S[co 2 ,li] would bring results closer 
in line with the model-based concept of ‘equilibrium’ fast-feedback 
sensitivity. The above-mentioned issues with aerosol influences mean 
that it would currently be best for estimates from palaeodata to report 
both S[ C02 ,li] an d 5[ C0 2 j li,ae]- 

Table 2 lists example estimates for S following the main potential 
permutations of the definition of S in our approach (for detailed break- 
downs, see Supplementary Information). The first example uses records 
of palaeodata since 800 kyr ago. The second example uses the same input 
data series 6 * , but focuses only on the LGM; the contrast between exam- 
ples one and two thus highlights state- dependence. The third example 
lists estimates for S [C o 2 ]> S[co 2 ,li] and S[co 2 ,li,vg] from a more model- 
based assessment for the mid-Pliocene (—3-3.3 Myr ago) 13 , with 
AT = 3.3 K relative to the present and AR CQ2 = 1.9 W m -2 due to C0 2 
increase from 280 to 400 parts per million by volume (p.p.m.v.; ref. 44). 
The broad range of S values found within each example illustrates that 
comparison across different definitions unrealistically widens the range 
of values reported, notably towards high values, because omission of 
‘forcing’ due to the action of any slow feedbacks will cause overestima- 
tion of S (see also Fig. 3). 

For a first-order estimate of the range of S from palaeodata that 
approximates compatibility with the centennial timescale ‘equilibrium’ 
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Figure 3 | Evaluation of results from Tables 1 and 2. y-Axis labels refer to 
numbered rows in these Tables, a, Data summary by table row. b, Probability 
assessment using normal distributions (shifted where relevant), c, Probability 
assessment using lognormal distributions. S[x\ refers to the climate sensitivity 
as defined in detail by the subscript X in Tables 1 and 2. For b and c, we assume a 
relative uncertainty of 25% for entries that lacked uncertainty estimates in the 
source studies. In a, rows from Table 2 are identified with either ‘plei’ or ‘plio’ to 
distinguish between the past 800 kyr and the Pliocene entries, respectively. The 
colour coding refers to broad geological intervals, as shown in the key. Boxes at 
right indicate which conditions were explicitly accounted for; that is, as 
‘forcings’ (in the CQ 2 /GHG column, filled squares indicate GHG and open 


squares C0 2 ). Circles (data points in a) show central values where reported, 
error bars represent uncertainties as outlined in the Tables, at the 1 o equivalent 
level. Arrow (case 21) indicates a value reported only as >0.8 KW” 1 m 2 . Black 
dashed lines in a show 68% probability limits for all estimates that account for at 
least ‘C02’ and ‘IT, based on thick dashed lines in b and c, taking whichever 
68% value offers the widest (more conservatively estimated) margin. In b and 
c, the solid black line indicates the mode value (maximum), and the thin dashed 
lines the 95% probability limits. All distributions in b and c are given as 
individual normalized frequencies (grey lines), and as mean normalized 
frequencies (red line). 
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Figure 4 | Equilibrium response of the global temperature as a function of 
C0 2 concentrations, based on three different approaches, a, This work, 
using data from the late Pleistocene of the past 800 kyr (ref. 6). b, Using data of 
the past 20 Myr (RW_1 1; ref. 61). c, Based on JH_12 (ref. 85) using similar data 
of the past 800 kyr as in a. d, Combination of all three approaches. Plotted areas 
include uncertainty estimates of one standard deviation. Because this work and 
JH_12 developed their approach only on Pleistocene data (climate being mainly 
colder than today), extrapolation of the impact of slow feedbacks to 2 X C0 2 is 


not meaningful (we show only extrapolation with fast feedbacks). RW_11 in 
contrast also includes warmer climates with C0 2 up to 450 p.p.m.v., so that the 
applicable range with slow feedbacks extends to 450 p.p.m.v. For future climate 
with 2 X C0 2 and a short time horizon (<100yr), only fast feedbacks are of 
interest (see d). Approaches partly disagree because of different assumptions. 
Uncertainties in this work (a) are estimated to be larger than they were in 
RW_11 (b) and JH_12 (c). For details of the equations and values used, see 
Supplementary Information. 


values of the IPCC-AR4 1 , values need to be used that account for ‘C0 2 ’ 
or ‘GHG’ as well as ‘LI’, and preferably also ‘AE’ and/or ‘VG’ (Tables 1, 2; 
Fig. 3). Such an assessment, excluding the case of row 21 in Table 1, 
yields a likely 1 (68%) probability range of 0.6-1.3 K W -1 m 2 , and a 
95% range of 0.3-1. 9 K W -1 m 2 (Fig. 3). These represent the widest 
margins out of two assessments, using either normal distributions 
with shifts when relevant (Fig. 3a), or lognormal distributions that in- 
herently allow asymmetry 2 (Fig. 3b). These assessments include uncer- 
tainties as outlined in the source studies, as well as any unaccounted-for 
dependence on different background climate states, but exclude poten- 
tial additional uncertainties highlighted in this study. Inclusion of 
ESS values (approximated by S[ C o 2 ]) would extend the upper limit 
beyond 3KW” 1 m 2 (Fig. 3a). Future work following a strict framework 
for reporting and comparison of palaeodata may refine the observed 
asymmetry. 

Finally, following our conceptual framework, we can make a projection 
of equilibrium temperature change over a range of C0 2 concentrations 
while considering either slow and fast (or only fast) feedbacks (Fig. 4; see 
Supplementary Information for details). Including the known uncer- 
tainties associated with palaeoclimate sensitivity calculations, and com- 
paring with two previous approaches 61,85 , we find overlap in the 68% 
probability envelopes that implies equilibrium warming of 3.1-3.7K 
for 2 X C0 2 (Fig. 4), equivalent to a fast feedback (Charney) climate 
sensitivity between 0.8 and 1.0KW” 1 m 2 . For longer, multi- centennial 
projections, some of the slow feedbacks (namely vegetation-albedo and 
aerosol feedbacks) may need further consideration. However, their impact 
is difficult to estimate from palaeodata, because uncertainties are large, 
and because responses during climates colder than present may differ 
from responses during future warming. 


We have employed a new framework of definitions for palaeoclimate 
sensitivity. This reveals how a broad selection of previously published 
estimates for the past 65 Myr agrees on a best general estimate of 

0. 6. 1.3 K W -1 m 2 , which agrees with IPCC-AR4 estimates for equilib- 
rium climate sensitivity 1 . Higher estimates than ours may suggest different 
climate sensitivities during particular periods, but a considerable portion of 
the higher values may simply reflect differences in the definitions of 
palaeoclimate sensitivity that were used. 
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CORRECTIONS & AMENDMENTS 


ERRATUM 

doi: 10.1038/naturel 1828 

Erratum: Making sense of 
palaeoclimate sensitivity 

PALAEOSENS Project Members 

Nature 491 , 683-69 1 (20 1 2); doi : 1 0. 1 038/nature 11574 

The surname of author V. Masson-Delmotte was misspelled as 
Masson-Demotte. In Table 1 of this Perspective, the ‘Source’ for 
‘Label in Fig. 3’ entry 13 should be ref. 61 (not ref. 65). In addition, 
the x-axis labels of Fig. 3c were misplaced. Figure 1 shows the cor- 
rected panel and the original Perspective has been corrected online. 
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Figure 1 | This is the corrected Fig. 3c. 
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